ar X iv : 0 70 7 . 05 02 v 1 [ m at h - ph ] 3 J ul 2 00 7 DEFLATED GMRES FOR SYSTEMS WITH MULTIPLE SHIFTS AND MULTIPLE RIGHT - HAND SIDES
نویسنده
چکیده
We consider solution of multiply shifted systems of nonsymmetric linear equations, possibly also with multiple right-hand sides. First, for a single right-hand side, the matrix is shifted by several multiples of the identity. Such problems arise in a number of applications, including lattice quantum chromodynamics where the matrices are complex and non-Hermitian. Some Krylov iterative methods such as GMRES and BiCGStab have been used to solve multiply shifted systems for about the cost of solving just one system. Restarted GMRES can be improved by deflating eigenvalues for matrices that have a few small eigenvalues. We show that a particular deflated method, GMRES-DR, can be applied to multiply shifted systems. In quantum chromodynamics, it is common to have multiple right-hand sides with multiple shifts for each right-hand side. We develop a method that efficiently solves the multiple right-hand sides by using a deflated version of GMRES and yet keeps costs for all of the multiply shifted systems close to those for one shift. An example is given showing this can be extremely effective with a quantum chromodynamics matrix.
منابع مشابه
Deflated Gmres for Systems with Multiple Shifts and Multiple Right-hand Sides∗
We consider solution of multiply shifted systems of nonsymmetric linear equations, possibly also with multiple right-hand sides. First, for a single right-hand side, the matrix is shifted by several multiples of the identity. Such problems arise in a number of applications, including lattice quantum chromodynamics where the matrices are complex and non-Hermitian. Some Krylov iterative methods s...
متن کاملDeflated BiCGStab for linear equations in QCD problems
The large systems of complex linear equations that are generated in QCD problems often have multiple right-hand sides (for multiple sources) and multiple shifts (for multiple masses). Deflated GMRES methods have previously been developed for solving multiple right-hand sides. Eigenvectors are generated during solution of the first right-hand side and used to speed up convergence for the other r...
متن کاملGlobal GMRES with Deflated Restarting for Families of Shifted Linear Systems
Many problems in science and engineering field require the solution of shifted linear systems with multiple right hand sides and multiple shifts. To solve such systems efficiently, the implicitly restarted global GMRES algorithm is extended in this paper. However, the shift invariant property could no longer hold over the augmented global Krylov subspace due to adding the harmonic Ritz matrices...
متن کاملDeflation Methods in Fermion Inverters∗
I will review recent developments in matrix deflation methods, by Ronald Morgan/Walter Wilcox, Andreas Stathopoulos/Konstantinos Orginos, and Martin Lüscher, with application to lattice QCD fermion inversion. I will begin with a short review of deflation-related work in the field. The Morgan/Wilcox algorithms using GMRES and BiCGStab for deflation will be described. Typical results for quenched...
متن کاملar X iv : p hy si cs / 0 60 91 58 v 1 1 9 Se p 20 06 Channel plasmon - polaritons : modal shape , dispersion , and losses
We theoretically study channel plasmon-polaritons (CPPs) with a geometry similar to that in recent experiments at telecom wavelengths (Bozhevolnyi et al., Nature 440, 508 (2006)). The CPP modal shape, dispersion relation, and losses are simulated using the multiple multipole method and the finite difference time domain technique. It is shown that, with the increase of the wavelength, the fundam...
متن کامل